INFN-ISS 93/2 
WIS-93/50/Jun-PH 



Exact response of the non-relativistic 
harmonic oscillator^ 

E. Pace^, G. Salme^ 

^ Dipartimento di Fisica, Universita di Roma "Tor Vergata", INFN, Sezione Tor Vergata, 
Via della Ricerca Scientifica, 1-00133, Roma, Italia 
^INFN, Sezione Sanita, Viale Regina Elena 299, 1-00161 Roma, Italy 

A.S. Rinat 

Weizmann Institute of Science, Rehovot 76100, Israel 



Abstract: 

Using Green's function and operator techniques we give a closed expression for the re- 
sponse of a non-relativistic system interacting through confining, harmonic forces. The 
expression for the incoherent part permits rapid evaluation of coefficients in a l/g expansion. 
A comparison is made with standard approximation methods. 
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1. Introduction 

Consider a system with mass composed of A constituents with equal mass m. We 
are interested in its hnear response or structure function, when the system is probed at 
momentum and energy transfer {q,u!) 

Siq,u) = A-^ \^On{q)fS{oj - /2mA - Eno), (1.1) 

n 

where J^oniq) = {0\pl\n) is the inelastic form factor between the ground state |0) and excited 
states |n), and pq = X^^- e^e*^'^ the intrinsic charge density. Eno are intrinsic excitation ener- 
gies without target recoil, and Pj = fj - R intrinsic coordinates with R the centre of mass 
coordinate. 

A formal summation over n yields an alternative expression 

Siq, oj) = (nA)-^ ^ Im ({0\e-'^<G{uj + Eq - q'^/2mA - ir?)e'«-^i |0)) 

i,j=i..A (1.2) 

G{z) = (z-K-V)-'^ 

with G, the exact Green's function of the system in terms of the total kinetic and potential 
energy 

Over the years, general methods have been devised to approximately calculate the above 
response. Recently there has been interest in the high-5 response for systems with singular 
attractive interactions of the confining type^'^. For those systems, perturbative methods in 
the confining interaction V fail and non-perturbative methods are required. We first mention 
the 1/q expansion of the (reduced) response as discussed by Gersch, Rodriguez and Smith 
(GRS)3 

Hq, y) = {q/m)S{q, u) = Y^im/lf Fkiv), (1.3) 

k 

where the energy loss w has been replaced by a suitable alternative variable 

y=-q/2 + mw/q (1.4) 
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Alternatively one may perform an actual summation of the series (1.1) ^. A pre-requisite 
for such a treatment is the knowledge of all the states of the complete spectrum, which is 
purely discrete for confining interactions. 

None of the above methods leads to a closed expression for S or </). However, Eq. (1.2) 
shows that such an expression would result if the exact Green's function there were known. 
This is only the case for a few selected systems. An enumeration of cases, which can be 
solved by means of path integral methods have recently been given At this point we 
recall that the above method has been previously suggested and applied to the response of 
systems, as complicated as liquid ^He probed at high q for fixed y ^. Evaluations in practice 
will always require approximations. 

In what follows, we will limit ourselves to a system of particles interacting through 
harmonic confining forces, which can be exactly solved and for which approximations can 
be tested. Since for this system the Hamiltonian separates in Jacobi variables, it suffices to 
treat only a 'di-quark' and actually for the same reason, one in one dimension only. In the 
following, the mass ruA will be replaced by 2m. 

In Section 2 we will discuss two methods which lead to a closed expression for the reduced 
response for the above system. In Section 3 we will give its l/q expansion and compare results 
with the outcome of the GRS theory as well as with previously suggested approximations. 

2. The response for a harmonically confined 'di-quark'. 

2a. The Green's function method. 

Consider the Hamiltonian of the relative motion of a harmonically bound 'di-quark' 



with (3, the inverse harmonic oscillator length. For instance the path-integral method leads 
to the following expression for the Greens' function for the relative motion ^. 



H = f/m + (3^x^/m, 



(2.1) 




(2.2) 
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where a = ^fp-sjq ; s = tqjm. The ground state wave function is (a;|0) = {0^ j-K^l^e'^'^'^'^ 1"^ 
and the part of the charge density fluctuation dependent on the relative coordinate pq{x) = 
_|_ g2g-i9x/2 Replacing the relative coordinates x,x' by 

x = Z + z/2, x' = Z-z/2, (2.3) 

one finds 

p\{x)pq{x') = e^e-'«^/2 + e^e'«^/2 + 6162 (e'«^ + e-'«^) (2.4) 

The first two terms are (for equal charges identical) incoherent contributions due to a trans- 
fered q, which is absorbed and emitted by the same particle. For the remaining coherent 
contributions those particles are different. The same transformation (2.3) (which has nothing 
to do with a CM transformation, and is typical for the harmonic oscillator case) permits an 
easy evaluation of the Z, z integrals. Taking the Fourier transform of (2.2) for the argument 
u + - q^ /Am, as required in Eq. (1.2) and using (1.4) to eliminate lo in favour of y, one 
shows 



incoh 



(e? + ei) 



/32s2 ^ sina/2 
a/2 



exp 



-iqs/A 



sma 



a 



(2.5a) 



0(9,2/)''°'' = 6162 / ^exp[ii/s]exp 



0^s^ fcosa/2 



exp 



. ,^ , sma 
iqs/4 ( -— + 1 



(2.56) 



4 V a/2 

Except for the harmonic oscillator, we do not know of a system for which the response 
has been given in a simple closed form hke (2.5). As is the case with the GRS theory, the 
incoherent part (2.5a) is seen to permit a series expansion in a oc However, the GRS 
expression ^ comprises three integrals and in addition an ordering operator which can only 
be realized in the expansion (1.3). 

4 



2b. Derivation using operator techniques 

In what follows we will sketch a different approach to obtain the response S{q,u), Eq. 
(1.1), for the case of a 'di-quark'. We do so, by directly evaluating the time-depending 
operators after using the Fourier transform of the S function. Eq. (1.1) then becomes 



S{q,uj) = \j (0|pj(*) ^^V>{-iHt) pq{x) exp 



i{ui 



2mA 



Consider first the time-translated space coordinate x{t) i.e. 



x{t) 



i=l 



t=0 



where 



H)t 



|0) 



(2.6) 



(2.7) 



{-iT[H, [H,[H, [H, 



(2.8) 



n times H 

Using elementary commutation rules and the hamiltonian, Eq. (2.1), for the relative 
motion, one obtains 



x{t) = X cos{2l3^ — ) - ^sin(2/32 A) 



(2.9) 



We now apply the Glauber formula e^+^ = e^e^e~^[^'^l, which holds if [A,B] is a c 
number'^, to the operator exp{iqx{t) /2) 



exp(ig^) = exp|i| 



exp 



iq^cos{2(3^-] 
2 m 



2 

exp 



X cos(2/32 A) _ P sin(2/?2 — ) 
m m 



zg-^sin(2/32l) 
2/5"= m 



exp 



-z^cos(2/32-)sin(2/?2l) 



(2.10) 



Inserting in Eq.(2.6) the unity operator / dk/{2'K)\k){k\ = I, and replacing by 2m, one 
obtains 



S{uj,q) 



dt 
2^ 



exp 



- —)t 
Am 



/dk 
— (0|pj(x) exp(-iiJt) pq{x) eMiHt) \k} {k\0) 



-I 



dt 
2^ 



exp 



i{u: - - — )t 
^ Am 



exp 



'1" ..:^tAQ2 t 



16/3= 



:sin(4/?^-) 



m 



dk 
27r 



M±^(0|fc+|[co.(2/3^^)-l]) + 



+ eie2 (0|/c+|[cos(2/32l) + l]) 
2 m 



exp 



m 



(fc|0) 



(2.11) 



After Fourier transforming the ground state wave function one can easily perform the 
integration over k. Replacing the variables t and w by s and y, respectively, one obtains Eqs. 
(2.5a) and (2.5b). 

The approach of this section can be extended to any hamiltonian for which explicit 
expressions for (2.7) and the initial state are available. Examples are the harmonic oscillator 
in any excited initial state and many body-systems interacting through harmonic oscillator 
forces. 

3. 1/q expansion, standard approximations, conclusions. 

Let us expand the integrand of the incoherent part, Eq. (2.5a), as a series in 1/q. In 
order to enable the integration term by term, Eq. (1.1) (and therefore Eq. (2.5)) has to be 
convoluted with some suitable smearing function, for instance the experimental resolution. 
This procedure results in a finite asymptotic limit for g — > oo. In contradistinction the 
coherent part can be shown to vanish exponentially in in that limit and is totally negligible 
for finite, large q. 

One finds for the coefficient functions in the GRS expansion (1.3) 



Fo{y) = 



{m/q)F,{y) = - ('^^ (^l 



1 /2/3 



3 \(3 



My) 



((m/g)^F2(2/) = --(^^ 



y \ I y 



3 V/3 



Fo{y) 



(3.1) 



My) 



We now compare this result from (2.5) with expressions and approximations which have 
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routinely been used in the past. First comes the GRS expansion (1.3) ^. One checks that 
the first two coefiicients above agree with those given before ^. The expression for F3 is new 
and higher order coefficients Fk may be calculated as well. The following comments are in 
order : 

i) {m/q)^ Fk{y)/Fo{y) is a polynomial of order 3A: in 

ii) the relevant expansion parameter is 2/3/g, 

iii) For increasing q at fixed y the response acquires non-zero contributions from ever increas- 
ing values of the energy with n — » 00 for g — > 00. Since for n ^ 00 the harmonic oscillator 
wave function tends to a plane wave^, one obtains the standard asymptotic limit Fq. 

Next we mention the first cumulant approximation ^° and an iteration of Fi The 
former makes sense only for at least three constituents and we therefore discuss only the Fi 
iterant. With 0o(a;) = (a;|0), the ground state wave function 

(j){q,y) = (27r)"^ / dse'^" / dx^x - s)0o(a:)exp \{-im/q) f da[Vix -a)- Vix) 

J-00 J K Jo 

= (27r)-iy d,e^J/«exp[-i/32s2]exp[l^]exp[--^] (3.2) 

= (27r)-i r dse^y^exp ["l/^'s'l U + h \^] + ^(e"') j 
J— 00 Jl "L^J J 

Note that the coefficient of the term 1/q is equal to the one obtained from Eq. (2.5a), whereas 
the other ones are quite different. 

In conclusion we have applied the Green's function and an operator technique to obtain 
the response of systems interacting through harmonic confining forces. For those systems, 
one can give a closed expression, Eq. (2.5), which can be tested against approximations^'^°'^^. 

The method can also applied to other systems for which the Green's function in closed 
form is available, as is the case for a ^-shell potential and a particle in an infinitely deep 
square well^ which has been previously investigated, using the summation method^ 

Finally let us note that the form (1.1) shows the response to be a sum of 5-functions in 
the energy transfer w. Our approach starts with the everywhere regular Green's function in 
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the conjugate variable t. After performing the manipulations and integrations as prescribed 
by (1.2), one obtains the reduced response (i){q,y)i Eq.(2.5), which is the Fourier transform 
of a regular function of q and t. On the one hand all of the above equations are only valid in 
the realm of the distribution function theory. On the other hand the expansion (1.3) for the 
reduced response with the coefficients (3.1) is a regular function of q and y. This apparently 
occurs when w is replaced by y. the spacing between the roots of the arguments of the delta 
functions appearing in Eq.(l.l) vanishes for q^ oo. 
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